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This paper describes highlights of an ongoing validation effort conducted to assess the 
viability of applying a set of analytic point source transient free molecule equations to model 
behavior ranging from molecular effusion to rocket plumes. The validation effort includes 
encouraging comparisons to both steady and transient studies involving experimental data 
and direct simulation Monte Carlo results. Finally, this model is applied to describe features 
of two exotic transient scenarios involving NASA Goddard Space Flight Center satellite 
programs. 


Nomenclature 


Variables 



A\ 

= 

normalization factor 

D 

= 

generic exponent 

d 

= 

diameter 

f 

= 

molecular distribution function, 



also a correlation parameter 

g 

= 

gravitational acceleration vector 

M 

= 

Mach number, also a generic flux 

m 

= 

mass 

m 

= 

mass rate 

n 

= 

number density 

h 

= 

vector surface normal 

P 

= 

pressure, momentum flux 

Qi 

= 

point source description 

4 

= 

energy flux 

R 

= 

gas constant 

r 

= 

distance magnitude 

s 

= 

speed ratio pu e 

T 

= 

temperature 

t 

= 

time variable 

u e 

= 

vector bulk exit velocity 

Kv 

= 

vector velocity variables 

V 

= 

velocity magnitude 

w 

= 

SCOS0 

X 

= 

vector position variable 

y 

— 

SCOSfZig 

z 


a- w 

a 

= 

prlt 

P 

= 
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= 

specific heat ratio 

J(x) 

= 

Dirac delta function 

0 

= 

angle between u e &x 

P 

= 

density 

T 

= 

shear stress, also time constant 

<f> 

= 

angle between x & n 

<t>e 

= 

angle between u,, & h 

6 

= 

mass flux 

¥ 

= 

angle between v & target surface normal 

Subscripts 


AV 

= 

a reference value 

e 

= 

exit condition 

i 

= 

incident 

INT 

= 

internal energy component 

n 

= 

number density 

L 


laser irradiation 

orb 

= 

orbital condition 

r 

= 

reflected 

TR 

= 

translational energy component 

0 

= 

stagnation condition 

1 

= 

associative marker 

1 

= 

perpendicular to surface 

D 

= 

lunar condition 

oo 

= 

ambient condition 



I. Introduction 


T he analysis and simulation of gases expanding from sources such as rocket nozzles into vacuum, or the effects 
plumes from these sources create when they interact with solid surfaces, present a considerable challenge to the 
scientific and engineering communities. As a plume expands into vacuum, density levels, and hence collision rates, 
decrease rapidly by many orders of magnitude. The main difficulty lies in accurately describing a flowfield 
extending from continuum flow at the nozzle exit, through the transition regime, and reaching free molecule 
behavior within a relatively short distance downstream. 

For thrusters, flow at the nozzle exit is usually characterized by high exit velocities and relatively high Mach 
numbers. Even in regions where significant intermolecular collision rates occur, relative velocity levels are low and 
little thermal scattering occurs normal to the mainly radial streamlines. Such observations lead one to consider 
describing the expansion under certain circumstances using free molecule (FM) theory. 

Under such assumptions, a model was created by solving the Boltzmann equation for a point source with a 
Lambertian thermal velocity distribution superimposed on a bulk, convective exit velocity. This paper describes a 
number of comparisons of model results with experiment and direct simulation Monte Carlo applications, both 
steady and transient. Example cases involving thrusters, molecular beams, and pulsed laser ablation will be 
presented. Along with ongoing validation efforts, novel applications involving satellite programs at NASA Goddard 
Space Flight Center will be discussed. 


II. Model Development 

To develop this model, the collisionless Boltzmann equation had been solved for distribution function / due to a 
directionally-constrained point source Q\ meant to describe directed flow from a nozzle exit over 2 n steradians 
centered on the source normal. 1 
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Qt represents directed flow from a Lambertian, thermal velocity distribution superimposed on convective exit 
velocity « e for mass flow rate rh , P = l/ V 2 RT , and speed ratio s = P u e . The influence of gravitational 
acceleration g is typically neglected for the applications of interest based on characteristic length scales and time-of- 
flight estimates. Normal n represents the orientation of a local starting surface element, and v • it emphasizes the 
imposed directional constraint. Generally u c is not aligned with h , the angle between the two being defined by (j) c 
(Fig. 1). Local angle (f> is measured between variable position x (distance r, experiencing local velocity v) and ft , 
and angle 6 is measured between w c and x. 




target normal 


L orb 


. target 
location 


FIGURE 1. Schematic representation of various quantities and angles used in analytical model. 

A. Steady Source 

The distribution of various flowfield quantities is found by integrating successive velocity moments of the 
resulting distribution function. In particular, the steady-state density p, mass flux O , incident normal momentum 
flux p j_ , and translational energy flux <jr TR are given by: 
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The normalization factor may be written as 



■^y(l + erfy). 


(6) 


Eq. (4) becomes incident normal momentum flux p ± or shear stress r when multiplied by COS!// or sin t// , 
respectively. When the source-sink geometry is fixed, V orb = 0. In these equations, r and </> are components of 
position x, measured from the centerline at the nozzle exit (Fig. Y), z = a - w, a = P r/t , w = scos^, and 
y = s cos (jy . Gas properties u c , exit temperature T, and mass flow rate m are defined locally on the starting 
surface. For diatomic or polyatomic molecules with specific heat ratio y, the internal energy flux <y INT is given by: 


f?INT 


(*) = 


r 5-3y ^ 

r - 1 


cjj(jc 
4 p 2 


(7) 


In addition, Eqns. (2) - (4) may be combined to obtain local expressions for velocity v and translational 
temperature T’tr : 


v(x) = 


p{x) ’ 


T’tr (*) 


i 

3 R | p{x) 



(8) 

(9) 


The total plume gas heat flux incident on a surface is given by q t = </ TR +<7int • Reflected heat flux q r is 
assumed to be carried away by diffusely reflected or desorbed plume molecules. The density of material is found by 
assuming the steady flux of molecules leaving a surface is equal to their rate of arrival, 2 given by Eq. (3). One will 
recognize the resulting density distribution is consistent with Lambertian outgassing. After some algebra, one finds: 
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( 10 ) 


For reflected quantities, surface temperature T r is a transient quantity. Over time, the net heat flux should 
decrease as T r increases based on conduction and radiation heat transfer mechanisms. For operations having short 
duration, it may be reasonable to hold T r constant in order to compute pressure p using Eq. (4), noting p = p L + p r . 

Because FM models ignore intermolecular collisions, collisional effects will be absent from their solutions. 
These effects include ambient and self-scattered return flux, evolution of chemistry, and externally-developed 
thruster interaction shock waves. The flowfield emanating from this point source is only valid where x ■ h > 0 due 
to the velocity constraint in Eq. (1). This last limitation is overcome by using local conditions from a starting 
surface that captures the expansion of gases expanding around the nozzle lip. 

B. Pulsed Source 

When m in Eq. (1) is described by a Dirac Delta function A m§(t), one may solve for successive moments off 
to find pulse-mode response expressions for various fluxes of interest. 


( \ 2A»ia 4 cos(zi w 2 -p ^ z 2 
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pr 
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III. Model Development 

This model was initially used to help study Space Shuttle/Mir thruster interaction problems during docking 
maneuvers, where it provided checks on direct simulation Monte Carlo (DSMC) calculations. 3 More recently, it has 
been used in support of a number of NASA-GSFC and -LaRC projects, with many successful comparisons to 
DSMC and experimental data. 4 " 6 Applications include effusion (Mach number M = 0), sonic orifices (M = 1), and 
thruster plumes (M > 4). 

A. Sonic Orifice 

An experimental investigation had been conducted at NASA Goddard Space Flight Center (GSFC) to observe 
the angular distribution of molecular nitrogen venting through circular tubes of different diameters and mass flow 
rates. 47 The test was conducted at T 0 = 290 K using GSFC’s Thermal Vacuum Chamber 238 (C238), with 
background pressure levels /> C2 38 ' n the 10 -10" 3 torr range during operations and driving pressures varying from 
2-50 torr. These conditions ensured that the N 2 issuing from the end of the tube would achieve sonic conditions, and 
for the lower driving pressure cases, Pc 23 H should not have interfered with the sonic expansion. 

Measurements were made using a Baratron capacitance manometer fitted to flexible tubing, the open end of 
which was attached to a goniometer/crankwheel mechanism. It appeared that for cases having the lowest driving 
pressures ( »i N below 0.02 g/s, Pc23& < 4 x 10" 4 torr), data collected from the smallest cross-section duct (d = 0.48 

cm) might provide a suitable set of conditions for making comparisons with the free molecule plume model. Data 
for these cases meeting these criteria were normalized by their centerline values and used to develop a best fit 
angular distribution curve by the study’s author. 7 This fit had the following form: 

p occos 3 0. (16) 

A comparison of this angular distribution with plume model results for a single point source of molecular 
nitrogen at AT = 1 is presented in Fig. 2. 



NORMAL TO SURFACE 


FIGURE 2. Comparison of Ref. 7 thruster model to FM sonic orifice flow for molecular nitrogen, Eq. (4). 
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B. CONTOUR Backflow 

On 15 August 2002, the NASA/Johns Hopkins University Applied Physics Laboratory (APL) Comet Nucleus 
Tour (CONTOUR) spacecraft was destroyed during conclusion of a solid rocket motor orbit transfer bum. 8 A NASA 
Mishap Investigation Board (M1B) was established one week after loss of contact, and under its direction, a number 
of failure modes were investigated and revisited in greater detail. 8 Ultimately, heating of the spacecraft due to 
plume radiation and convective impingement in the backflow regime was reasoned to be the first element of the 
proximate cause for its demise. 8 In the course of the investigation, a computational fluid dynamics (CFD) solution 
for the near-field gaseous expansion produced by NASA-Marshall Space Flight Center (MSFC) was used to provide 
a starting surface for DSMC and free molecule plume modeling efforts to determine fluxes in the backflow regime. 9 

A comparison of model flowfield results is presented in Figs. 3a & 3b. In the figures, the CONTOUR spacecraft 
(S/C) forebody is shown on the left, while the CFD starting surface is represented by the shadowed region to the 
right. In between are flowfield density contours for the plume impinging on the satellite from right to left. These 
exhibit a high level of agreement over three orders of magnitude. The free molecule contours consist not only of the 
direct contribution from the plume (at high speed ratios s), but surface reflected molecules as well ( s = 0). 



FIGURE 3. Comparison of steady DSMC/DS2V (a) and FM results (b), 
logarithmic number density contour maps. 


Sampling the flowfield axially at different radii revealed the high levels of agreement were maintained nearly 
everywhere, even by the S/C Payload Attachment Fitting (SPAF), except very close to the spacecraft forebody 
surface. In the DSMC solution, a radial Knudsen layer had developed where density increased as the surface was 
approached from a normal direction. In the free molecule model, mass conservation ensured that density levels 
would drop instead of increase there. Incident heat flux profiles across most of the forebody were similar, except the 
presence of a SPAF attachment ring in the DSMC geometry led to a localized intensity of heating that was absent 
from the free molecule calculation. 


C. Molecular Beam Time-of-Flight Data 

One novel feature of this plume model is its transient flow capability. When considering model validation, it is 
difficult to recreate highly-rarefied conditions using ground-based facilities due to demands for sufficient pumping 
capacity to eliminate the influence of vacuum chamber surface scattering on a thruster plume expansion. 1011 
However, it is possible to compare limiting cases of such models to data collected during carefully executed 
experiments employing molecular beams (MB). In particular, a number of early experiments measured molecular 
time-of-flight (TOF) data to produce velocity distributions associated with low-density sources. 12 13 

In one such experiment. Miller & Kusch measured TOF data for effusion (M = 0) of potassium and thallium at 
900 K and confirmed the on-axis distribution of density was given by: 
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( 17 ) 


Such a relationship is also given by Eq. (11) for pulse-mode flow for speed ratio s = 0 and t/>=0 (Ref. 1 1). 

In another molecular beam experiment incorporating a nozzle for their source, Anderson & Fenn found for high 
Mach number TOF measurements along the nozzle centerline: 13 

2 

pxa 4 e~ = . (18) 

Again, this relationship is consistent with Eq. (1 1), but for s > 0. 

D. Transient Pulsed Laser Ablation Simulation 

Bykov et al. modeled a transient pulsed laser ablation process for heating and vaporization of a niobium target 
with varying laser fluence levels (4-6 J/cm 2 ) of brief duration, using direct simulation Monte Carlo (DSMC) for the 
subsequent vapor release. 11 14 It was noted that the response of vaporized effluent angular distributions to increasing 
heat inputs was to change from a spherical expansion to one more focused along the surface normal initially. Such 
distributions would relax to spherical expansions as time and rarefaction levels increased. 14 A plot of centerline 
density function f n versus a was presented (Fig. 4), where 


fn - 


n t 


3 


n AV r L 


(19) 


In Eq. (19) n is the local number density, n A y is a reference value based on the total amount of removed mass 
(dependent upon laser fluence level), and is the 13 ns laser irradiation period. 14 The f n curve for each laser 

fluence level features a maximum ( a ~ 0.8 for 4 J/cm 2 , 1.15 for 6 J/cm 2 ), and for large levels of a (short time 
periods) the functions decay rapidly. 



FIGURE 4. Comparison of Eq. (11) development with Ref. 14 DSMC simulations of transient PLA plume 

centerline responses for different heat flux inputs. 

From Eq. (11), one can write an equivalent centerline variation: 


nr oc f n oc a e 


( 20 ) 


which is maximized for 
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( 21 ) 


s + "\/s 2 +2 

a = 

2 

Not enough information is given in Ref. 14 to complete a direct comparison to their f n curves, but Eq. (21) 

indicates maxima shift from a= 1 / -Jl for effusive flow to larger values for s > 0, which would be consistent with 
the sort of energy transfer into more directed flows Ref. 14 authors observed for higher laser fluence levels. 
Substituting s = 0.33 & 0.715 in Eq. (20), corresponding to a values similar to those producing maximum f n values 
in the DSMC results listed above, model results provide fair matches for the 4 & 6 J/cm 2 cases with arbitrary density 
scaling. This level of agreement extends over a wide range of time scales (<x 1 /a) even though intermolecular 
collisions are responsible for deviations from effusive flow close to the model surface. 

IV. Other Applications 

Over time, ongoing validation activities such as the examples included in the previous section have led to 
increasing confidence in using the molecular flow plume model as a tool to describe an unusual variety of plume 
interaction issues, particularly for NASA Goddard Space Flight Center spacecraft missions. Two applications 
featuring transient effects are briefly presented below. 

A. CloudSat-CALIPSO Interactions 

CloudSat and CALIPSO are the second and third elements of the Afternoon Constellation, or “A-Train”, a joint 
effort between NASA and ESA to operate a heterogeneous group of scientific satellites in a coordinated effort to 
make a variety of near-simultaneous atmospheric measurements over an identical ground track. 15 16 The A-Train 
operates in a sun-synchronous orbit at 705 km altitude. Each satellite’s Attitude Control System (ACS) features 
either 1 lb f or 1 N aft- facing monopropellant hydrazine thrusters. Thruster firing periods range from 7 - 100 s. 
Plume species were based on 1 00 percent decomposition of high-purity grade monopropellant hydrazine with 65 
percent ammonia decomposition, 17 adding 0.5 percent mass fraction of water vapor as a system impurity. A nozzle 
exit velocity V e = 2.2 km/s was assumed, with a mass flow rate m of 2.0 g/s per lbf thrust. 

TABLE 1. CloudSat Model ACS Plume Species Mass Fractions. 


Species 

Mass 

Fraction 

NH, 

0.247 

n 2 

0.667 

h 2 

0.081 

h 2 o 

0.005 


The MSISE-90 atmospheric model 18 generated estimates of atmospheric density for A-Train conditions of n K ~ 
2.1xlO u molecules/m 3 (±5x), essentially made up of atomic oxygen. The calculated mean free path length 
associated with plume-ambient molecular scattering under such conditions is immense (-7000 km). This estimate is 
only representative at an altitude of 705 km; clearly there will be an overwhelming level of interaction in the lower 
reaches of the atmosphere, and the plume will not develop symmetrically with altitude. 

Nevertheless, some program managers were concerned about the possibility that thruster operations by one 
satellite in this formation may create unwanted impacts on the function of others. CALIPSO trails CloudSat 
nominally by 15 s, equating to a separation distance of roughly 112.5 km. CloudSat employs a pair of 1 N thrusters, 
and ACS operations typically last 19 s. Due to the separation distance, peak impingement fluxes on CALIPSO 
should be experienced roughly 50-70 s after a CloudSat ACS unit commences firing. 

Due to the short ACS firing intervals relative to the time taken to travel distances separating A-Train spacecraft, 
it might not be adequate to assume steady conditions for determining surface flux couplings. It was decided to use 
Eqns. (2) - (10) in transient mode as a square wave to observe the development of CALIPSO impingement fluxes 
due to CloudSat ACS thruster operations. In addition to neglecting scattering occurring via ambient conditions at 
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705 km altitude, the effects of atmospheric density stratification with altitude and solar ionization have been 
neglected, making this study a first approximation to actual impingement behavior. 

IE-12 
IE-14 
IE-16 
IE-18 

M | IE-20 

g ,E ' 22 

£ IE-24 

<z 

IE-26 

2 

IE-28 
IE-30 
IE-32 
IE-34 

1 10 100 1000 
TIME ELAPSED FROM INITIATION [s] 

FIGURE 5. Transient mass flux experienced by CALIPSO due to single CloudSat ACS operation. 

Results for total impingement mass flux ® (combination of N 2 , H 2 , NH 2 , and H 2 0 as described in Ref. 17) 
experienced by CALIPSO due to CloudSat operations are presented in Fig. 5 above. The initial steep rise in ® is 
due to dominance of the directionally-independent wave-like term a = fir/t in parameter z. However, as time 
progresses, the effect of bulk motion being directed downstream in the model causes a shift in this rise, manifested 
by the kink occurring at about 13 s elapsed time. 

Another deviation in the transient mass flux impingement profile occurs around the period anticipated for peak 
fluxes to occur. This is due to a species separation effect. The transient effective molecular weight of the gases 
reaching CALIPSO ’s ram- facing surface is presented in Fig. 6 below. 




TIME ELAPSED FROM INITIATION [s] 


FIGURE 6. Transient average molecular weight experienced by CALIPSO due to single CloudSat ACS 

operation. 

Before the pronounced bulge in Fig. 5 occurs, the gas consists almost entirely of molecular hydrogen. At high 
degrees of rarefaction, lightweight species will tend to spread fairly diffusely while heavier species remain more 
tightly focused around the plume axis, traveling more closely to the average bulk velocity. 19 The peak effective 
molecular weight reaches 22.7 g/mole at 64 s elapsed time, versus an expected average molecular weight of 12.9 
g/mole. At later times, heavier gases disappear and the residual ensemble consists of molecular hydrogen again. 
Notice how this gas appears to linger in Fig. 5 even though the thruster was only on for 19 s. At its peak, the effect 
of firing two CloudSat thrusters over a 19 s interval produced a density of 7 x 10" 19 g/cm 3 , an increase of roughly 12 
percent over ambient levels. 

Other observations for this application involved comparisons between mass flux and other fluxes. 15 Among 
other things, higher velocity moment solutions featured steeper final decay rates as square wave solutions 
degenerated to those for a pulse source, Eqns. (11) - (15). In the pulse source solutions, every increase in velocity 
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moment includes another power of f\ It was also noted that higher velocity moment peak fluxes occurred at 
slightly shorter elapsed intervals than lower ones. 15 

B. LADEE - Surface-Scattered Flux Environment 



FIGURE 7. Artist’s rendering of the Lunar Atmosphere Dust Environment Explorer. 

The Lunar Atmosphere Dust Environment Explorer (LADEE) spacecraft is being designed for a mission 
featuring low altitude orbits of the Moon to take relevant ambient measurements before that environment becomes 
altered by future activities associated with human exploration (Fig. 7). Depending on orbital insertion timing, the 
observatory’s mission could last over four months, beginning in 20 13. 211 In some cases, the orbital altitude may dip 
as low as 20 km above the lunar surface. 20 LADEE’s instrument suite includes a sensitive Neutral Mass 
Spectrometer (NMS) which will be oriented in the spacecraft’s ram direction to measure tenuous atmospheric 
constituents. Ambient density levels vary widely based on local solar illumination, but may exceed 10 5 
molecules/cm 3 during a lunar dawn. 21 NMS is being designed to detect individual species concentrations at much 
lower levels, however. 

LADEE’s attitude control system features 5 -lb f bipropellant thrusters mounted in two canted pairs beneath the 
spacecraft’s lower deck. 20 If the ACS units are directed forward and toward the lunar surface, it is possible for 
individual species associated with their operations, reflected from the daytime lunar surface, to be quantitatively 
measured by the NMS as the observatory subsequently passes overhead. At lower altitudes, thruster plume fluxes to 
the lunar surface are maximized, as are reflected fluxes back to the spacecraft; also NMS may intercept peak species 
fluxes before lunar gravitational disturbances exert much influence. In this manner, NMS measurements could be 
used to test a variety of characteristics associated with this transient plume model, including evolution of species 
concentration (magnitudes, time-of-flight measurements, angular distribution, species separation effects), thruster 
exit conditions, or to infer local elements of lunar surface properties such as permeation or gas-surface interaction 
models. 

LADEE will orbit the Moon in a variety of orientations and orbit altitudes. With surface values of 1737 km for 
lunar radius R and a gravitational acceleration constant of g = 1.62 m/s 2 , the orbital velocity at 20 km circular 
altitude should be F orb = 1.67 km/s. A study was developed" assuming a pair of these ACS thrusters facing 
forward and 20° below the satellite velocity vector were operated for just one second using the pulsed source model. 
Under these conditions, it was considered that gravitational effects on the gas at this altitude and for an ACS bulk 
velocity of F e « 3.0 km/s could be neglected. 22 Expected plume product nozzle exit temperature was T e « 550 K, 
and ACS plume product species are presented in Table 2 below. 2 ’ Interactions with a featureless lunar surface were 
modeled as impermeable, chemically-inert, and diffusely-reflecting. Such details allowed use of the plume model 
for s = 0 conditions and a simple mass conservation statement for the surface-reflected elements. 
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TABLE 2. Dominant ex pected LADEE ACS plume species mass fractions. 2 


Species 

Mass 

Fraction 

n 2 

0.426 

h 2 o 

0.294 

CO 

0.178 

o 

o 

0.086 

H 2 

0.016 


Under such conditions, results indicate LADEE’s NMS could indeed detect surface-reflected ACS thruster plume 
gases (Fig. 8) with peak values occurring within a minute after that brief operation. Individual species responses are 
separated by molecular weight, with peak molecular hydrogen measurements in particular occurring appreciably 
more quickly than expected based on bulk velocity and geometric configuration (19-20 s). Reviewing the pulse 
source equations and noting that flux M varies as t D . the time to reach a peak flux becomes 

^ max flux — 

W 1 

V 


2 ft r 


+ Jl + : 


2D 


■\ ' 


( 22 ) 



FIGURE 8. Predicted ACS species contributions measured by LADEE NMS, assumed conditions. 

For w = .y on the plume centerline, f maxflux — » rj V c as s —> oo . For any finite value of s however, the time to 
reach maximum flux actually decreases due to the component of kinetic energy associated with the non-zero thermal 
distribution. In Fig. 8, the peak surface-reflected H 2 concentration reaches NMS in less time than one would expect 
the average molecule to reach the surface! This effect was investigated further by artificially reducing T e to 55 K, 
which led to a tighter focus of gas tending to reach the lunar surface around 19 s and a delayed response reaching 
NMS compared to Fig. 8. 22 


V. Concluding Remarks 

A number of results have been presented from an ongoing validation effort for a transient free molecule thruster 
plume model over a wide range of applications. Encouraging results have led the author to assist NASA Goddard 
Space Flight Center and other organizations assess increasingly complex scenarios associated with exotic scientific 
spacecraft missions. It is hoped that such efforts will allow these organizations to perform their activities with 
greater confidence and possibly enhance general understanding. 
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Introduction 


Analysis and simulation of gases expanding into high 
vacuum present considerable challenges to engineering and 
scientific communities 

Examples for spacecraft environmental issues 

- plume impingement (forces, heating) 

- contamination of sensitive surfaces 

• optical instrument apertures (mass flux, density) 

• thermal control surfaces (mass flux) 

Flowfield passes from continuum to ffee-molecule regimes 
in relatively short distance 



Introduction (2 of 2) 


• Existence of collisions may not substantially alter plume 
distribution from free-molecule description 

- collision rates fall rapidly with distance from source 

- majority of self-scattering collisions occur when faster 
molecules overtake slower ones on similar trajectories 
from source 

- center-of-mass motion remains unchanged during 
collision 

• Can consider modeling this regime using free molecule flow 
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Plume Model Formulation-Source 


• Find particular solution to collisionless Boltzmann equation 
for source Q x : 


where 



and 






cos z q L 


+ s cos <j) e (l + erf (5 cos )) 
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Plume Model Development-Source 



• Impingement angle with target surface given by i// 

• Simplifies for axisymmetric conditions 

- 4 = 0 

- ( j )=0 

• other definitions: 

s = Bu P = , ; z = a-w ; a = Br!t ; w = scos<9 

J2RT , L 
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Model Development— Velocity Moments 


p{x,t) 


fim cos </> w 2 _ s 2 



+ 



+ w‘ 




J 


n erfc z 




rh cos 0 cosy/ 

A x 7ir 2 



\cc + CCW + W + 1^6 



(3 4 


+ 


4ft werfcz 


v 2 J 

J 


• Similar, but more complicated expressions for pj_,T, g TR > ^int 

• Can also estimate 

- diffuse, surface-reflected quantities ( s = 0) 

- local average velocity, <GMW>, temperature with suitable ratios 
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Model Development— Limitations 


• No inherent mechanism for gaseous collisional effects in free 
molecule flow 

- Molecular scattering 

- Chemical reactions 

• In present application 

- Neglect stratification of atmospheric density with altitude 

- Possible influence of solar ionization 

- Neglect physical presence of satellite bodies 

• Consider this study as first approximation 

- Could possibly be revisited if motivation existed 
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Model Development — Pulse 

• In plume model analytical development for this 

application, can also create velocity moments for case 
where mass flow rate is described by m — A mS{t) 

p(x,l)= 2Ama ‘‘ cos « > < ,-(»-») 2 < ,--- 2 . 

A\7ir 

pr^ 

0(x,f) = — cos^/; #tr {x,t) = —^cosy/; 

t/ ^ 1 

2 2 

( \ P r 2 / \ pr . 

Pi_\ x ->t) = cos y/\ T{x,t) = ^-smi//cosi// 

t t 
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Sonic Orifice Comparison 


• Test performed at NASA-GSFC to measure density 
angular distribution 

- N 2 at various p 0 , tube diameters 

• Appears conditions closest to FM flow for lowest p 0 , d { 

- Results suggest 

p oc cos 3 6 

• Can compare to steady plume model result for M = 1 
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Sonic Orifice Comparison 


COMPARISON OF ANGULAR DISTRIBUTIONS 
FOR MOLECULAR FLOW EXPANSIONS 
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CONTOUR SRM Backflow 

• CONTOUR spacecraft destroyed during conclusion of SRM 
orbit transfer bum; Mishap Investigation Board formed to 
determine proximate cause for this outcome 

• NASA-MSFC supplied CFD starting surface for DSMC 
study of high-angle satellite surface convective heat input 

- Simplified (gases only) 

- Asked to compare using plume model 

• FM flowfield result of superposition 

- Free expansion 

- Surface-reflected gases ( s = 0, no deposition) 
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CONTOUR Model SRM Flowfields 
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• Density contours, logarithmic scale 
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CONTOUR Model Comparison — Cmts. 


• Decent level of agreement over four orders of magnitude 
within flowfield 

• Comparisons of surface fluxes also exhibit good agreement 

- Pressure 

- Convective heat flux levels 

• Growth of Knudsen layer in DSMC model responsible for 
subtle differences 

- DSMC contours increase as surface is approached closely; 
FM flowfield contours fall to free expansion values 

- Buildups in surface fluxes displaced radially outward 
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Molecular Beam TOF 


• Transient model comparisons difficult to verify experimentally 

- Difficult to eliminate effects due to vacuum chamber 
surface scattering for thruster plumes 

• Exception: STG Facility, DLR-Gottingen 

• Can compare limiting cases for plume model with carefully 
executed molecular beam (MB) experiments 

- Time-of-flight (TOF) data used to produce velocity 
distributions associated with low-density sources 

- Miller & Kusch, effusion ( s = 0) of potassium and thallium 
@ 900 K 

- Anderson & Fenn, high-M nozzle source, (conditions?) 
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Miller & Kusch TOF Results 
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MB TOF Results 


• For effusion, Miller & Kusch determined 

4 -a 2 

pec a e 

• For high-Mach number, Anderson & Fenn found 

A 2 

4 — z 

pec a e 

• Both developments are reproduced by pulse source model 
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Pulsed Laser Ablation Simulation 


• Bykov et al. modeled pulsed laser ablation (PLA) process 
for heating and vaporization of niobium due to varying 
laser fluence levels (4-6 J/cm 2 ) 

- Used DSMC to simulate subsequent vapor release 

• Noted that effluent angular distributions changed from 
approximately spherical expansions to ones more focused 
along surface normal for early periods 

- Relaxed to spherical expansions as time and rarefaction 
levels increased 

• Presented centerline density using 



n AV T L 
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PLA - Model Development 


• Pulse mode model suggests that 



2 

ccf n ccae z 


• Maximized for 

s + a/s 2 + 2 
a = 

2 

- For effusive flow, occurs for a = 1 / V 2 

- For non-zero 5, critical a increases 
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PLA Model Comparison 



a 
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Outline 
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“Afternoon Constellation” 

• International effort to operate heterogeneous satellite group to 
make near-simultaneous measurements over defined ground track 

- Was to feature seven satellites, including OCO and Glory 

- Sun- synchronous orbit 

- 705 km altitude 

- One minute separation « 450 km 

- Each S/C features aft- facing N 2 H 4 Attitude Control System 
(ACS) thrusters 

• Atmosphere (MSISE-90 model) 

- 99+% atomic oxygen 

- n m « 2.1 x 10 11 /m 3 (± 5x) 
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A-Train Configuration 
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Satellite Formation Flying 
(FF) Interaction Example 


• Program managers sometimes voice concern that thruster operations 
from one spacecraft may create non-negligible effects on others 

- Schedule and cost always at a premium 

• Combination of thruster pulse durations and separation distances may 
mean impingement events should be modeled in a transient mode 

• In particular, estimate CloudSat ACS operation effects on CALIPSO 

- 1-N N 2 H 4 units operating for 19 s (step functions, square wave) 

- CALIPSO trails by 15 s (« 110 km) 


27-30 June 2011 


25 


MSW 



FF Thruster Parameters 


• A-Train S/C ACS feature sets of aft-facing monopropellant 
hydrazine thrusters 

- Mass flow rate = 2.0 g/s per lb f thrust, V e = 2.2 km/s 

• Plume species based on 65% NH 3 decomposition 

- H 2 0 included as fuel impurity at 0.5% mass fraction 

• Neglected 

- Gravitational effects 

- Scattering by ambient atmosphere 

• CALIPSO should encounter peak impingement fluxes 
roughly 50 - 70 s after CloudSat commences a 19 s operation 
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NORM. MOM. FLUX [Pa] DENSITY [g/cm 3 ] 


CloudSat/CALIPSO- 
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Various Fluxes 



10 100 1000 
TIME ELAPSED FROM INITIATION [s] 


• Notice post-peak slope 
increase with higher 
velocity moments 
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Example— CloudSat/CALIPSO, GMW 
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Induced Environment Effects 


• Depending on circumstances, thruster operation induced 
environments might not be insignificant 

- CloudSat — >CALIPSO (At = 15 s), single thruster: 

—19 3 

Anax = 3-5x10 g/cm 

• Two thrusters — >peak density ~12 percent of ambient levels 
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LADEE Introdction (1 of 2) 


• Lunar Atmosphere Dust Environment Explorer (LADEE) 

- Collect data regarding lunar atmosphere (gases, dust) before 
alteration due to future exploration activities 

• Features include 

- Operational period -100 days 

- Variety of orbits (elliptical, circular) 

• Nominal = 50 km, circular 

• As low as 20 km, circular 

- Variety of orientations used for making measurements, 
communicating with Earth 

• Lunar atmosphere is so rarefied it’s referred to as an “exosphere” 

- Essentially free-molecule conditions 

30 



LADEE Introduction (2 of 2) 


Instruments include Neutral Mass Spectrometer (NMS) 

- Design is sensitive enough to detect --100 molecules/cm 3 

NMS measurement sensitivity drives many LADEE 
contamination control requirements 

- Causes consideration of unusual scenarios 

• Outgassing 

• Attitude Control System (ACS) thruster plume influence 



LADEE Spacecraft 
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Schematic Diagram for Study 
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LADEE Case Study Objective 

• Appeared conceivable NMS could measure ACS plume gases 
reflected from sunlit, impermeable lunar surface 

- At minimum altitude 

• Measurement would be maximized 

• Gravitational influence minimized (“short” time-of-flight situation) 

- Could use to verify aspects of thruster plume modeling 

• Model the transient disturbance to NMS measurements due to 
ACS gases reflected from lunar surface 

• Observe evolution of various model characteristics as 
measured by NMS 

- Species magnitudes, TOF measurements, angular 
distribution, species separation effects 
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LADEE Test Case Conditions 


• Minimum altitude (20 km, circular) 

• NMS faces ram direction 

• Forward- facing ACS thruster pair 

- Operates for 1 s, 20° below horizontal 

- Ignore changes in spacecraft altitude (use pulse mode) 

• Units: 5-lb f MMH/MON-3 

- V e » 3.0 km/s, T e » 550 K 


Species 

n 2 

h 2 o 

CO 

co 2 

h 2 

Y, 

0.43 | 

0.29 

0.18 

0.086 

0.016 


• Particularly interested in water vapor influence 
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LADEE — Lunar Impingement: 

Surface Density 


20 s 


40 s 


60s 
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LADEE — Lunar Impingement: 
NMS Species Density Estimates 
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Nozzle Exit Temperature T e Influence 


Elapsed time for peak species mass fluxes to reach lunar surface 
occurred quicker than expected based on V e sin 20° 

Time derivative of mass flux equations (<D cc t~ D ) indicates 


t 


2f3r 


max flux 


f 


W 


1 + . 1 + 


2D 


\ 


V 


2 , 
W J 


For w = s on the plume centerline, t max flux — > r/ V e as s — > co 
For finite s, this period is always shorter 
- Consequence of thermal energy component 



Concluding Remarks 


• A FM model has been created for estimating thruster plume effects 

• Ongoing verification effort for a wide variety of conditions 
demonstrates utility of this approach 

• Results characterized by 

- Transient effects 

- Superposition 

- Species separation 

- Long-term persistence of lightweight species 
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Extra Materials 
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Example — Approach to Steady State (cont.) 


• Upstream locations tend to experience longer transients 
than downstream 

- Review variable z 

• Highest angles produced upstream, lowest cos lvalues, 
longest transient influences 


z=a-w- 



— s cos 0 


t 
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MASS FLUX [g/cm /s] 


Example — Impingement Fluxes 

ENHANCED FLYING FORMATION THRUSTER INTERACTIONS 
TOTAL MASS FLUX EXPERIENCED BY TARGET SPACECRAFT 



10 100 1000 10000 

TIME ELAPSED FROM INITIATION [s] 
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Example — Impingement Fluxes, <GMW> 

ENHANCED FLYING FORMATION THRUSTER INTERACTIONS 
AVERAGE MOLECULAR WEIGHT EXPERIENCED BY TARGET SPACECRAFT 



10 100 1000 10000 


TIME ELAPSED FROM INITIATION [s] 


27-30 June 2011 


43 


MSW 



Example — Impingement Fluxes (cont.) 


• Compare OCO & Aura environments 

- Early behavior dominated by wave-like term a 

• Essentially all H 2 

- Diverges due to effect of bulk motion downstream 

• Absent upstream 

• Still essentially H 2 

- Extra downstream hump due to increasing presence of 
other plume species 

• Maximum fluxes occur near time expected due to bulk velocity 

• Only H 2 left lingering at longest times 
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Example — Pulse, Step Comparison 


Click 

Contour 

Map 


Point 1 s Mass Flu* | 05 May2006 I GLORY A-TRAIN (ENDCAP) 


Time = 5 s 


1-lbf N2H4 AGS THRUSTER 
Time on = 1 s 
Point Source Soln. 
Mass Flux [g/cm2/sl 
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Example — Approach to Steady-State 


Click 

Contour 

Map 


Step.Steady I 05 May 2006 I G LORY A-TR AIN (E ND CAP ) 


Time = 5 s 
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>■ 
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Transient to Steady 
Step Source Soln. 

Mass Flux [g/cm2/s] 
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Example — 100 s Operation 


Click 

Contour 

Map 


10 0s Pulse | 05 May 20 06 I G LO RY A -TRAIN (E ND CAP ) 


Time = 5 s 


1.5E+08 - 
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>■ 
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Transient 100s Pulse 
Mass Flux [g/cm2/s] 
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K.E. FLUX [W/cm ] 


CloudSat/CALIPSO— Energy Fluxes 



TIME ELAPSED FROM INITIATION [s] 



TIME ELAPSED FROM INITIATION [s] 


• Long-term, model predicts eventually higher IE flux than KE flux 

- IE flux created as product of IE/mass and mass flux 

- Not actually realized (radiative loss of IE with time) 
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<V> [m/s] 


CloudSat/CALIPSO— Velocity 



• Not exactly given by r = vt\ deviations evident around V cxj 

- Effect of finite pulse length still significant 

- Superposition effect due to multiple species 
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